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Abstract 

We present a model for the formation of relativistic global vortices (strings) at 
a quench, and calculate their density and correlations. The significance of these 
results to early universe and condensed-matter physics is discussed. In particular, 
there is always open, or infinite, string. 



1 Introduction 



Many systems produce topological defects, in the form of vortices or monopoles, on 
undergoing a phase transition to an ordered state. In this paper we shall attempt to 
calculate how such topological debris is produced in simple transitions. 

As might have been anticipated, the distributions that we shall find are largely generic, 
roughly independent both of initial conditions and of the way in which the transition is 
implemented. It is for this reason that, hitherto, it has not been thought necessary to 
derive defect distributions in any detail. In fact, given our poor understanding of the 
underlying theories, it has been a consolation to be able to fall back on generic scaling 
solutions. For example, the large-scale structure of the universe has been attributed 
|H, to cosmic strings (vortices in the fields) formed at the Grand Unification era. 
Although we have only a primitive understanding of the relevant field theories, it has 
been suggested that the initial conditions of any string network are largely washed out 
after a few expansion times, at which the network is assumed to approach a scaling 
regime with a few large loops and long strings per horizon volume continuing to produce 
smaller loops by self and mutual intersection. 

Despite this, there are two circumstances in which this lack of detailed initial infor- 
mation leaves us at a loss, to which this paper is largely addressed. The first, of less 
interest, concerns the density of the defects formed. We have known for some time how 
to make reasonable qualitative estimates [|TJ of densities and in this paper we can confirm 
them quantitatively. Secondly, it is not entirely true to say that the details of the string 
distributions are unimportant. In particular, for the astrophysicists the scaling solutions 
for the early universe mentioned above arguably require some open or 'infinite 7 string, 
i.e., vortices that do not self- intersect. The presence of such string is, in part, determined 
by the initial conditions and it is important to know whether it is present for reasonable 
models and if so, how much. For example, it has been suggested that vortices produced 
by bubble nucleation in a strong first-order transition will only form small loops 0. 

The paper is organised as follows. We begin by reiterating the main tactics for deter- 
mining defect densities and distributions and then display the forms that they will take 
in a Gaussian approximation for the underlying field distributions. The main tools are 
the correlation functions of the defect densities. As it stands, some work is necessary to 
convert them into easily measurable, or easily identifiable, quantities. We present exam- 
ples, motivated both by our dynamical model (but simpler) and by current numerical 
simulations in astrophysics, to help us understand them better. 

In our model these correlation functions are given a concrete realisation in terms 
of exponentially growing unstable modes when the phase transition is implemented by 
an instantaneous quench. Although this is an unrealisable idealisation, these results 
are used as a benchmark for more general transitions when, later, we vary both the 
initial conditions and the way in which the quench is implemented. It was already 
implicit in our earlier work|4|], of which this is a continuation, that for a very fast quench 
the initial conditions will, in general, only determine the subleading behaviour of the 



1 



defect production. We shall extend the work presented there to look for exceptions to 
this general behaviour. Further, it will be seen that changing the rate at which the 
quench is implemented can be approximately equivalent to changing the time at which 
the transition can be said to have begun. As a result, the preliminary work of |3J for 
instantaneous quenches can be extended with only minor modification. 

After a discussion of the way in which defects freeze into the field we conclude with 
some observations about the length distribution of vortices, obtained by interpreting 
some recent numerical simulations |5|, |6| in the light of our model. It will be seen that, 
whatever we do, there will be infinite string. Yet again we assume flat spacetime, for 
simplicity. Our conclusions are thus, in this regard, more applicable to weak coupling 
condensed matter physics than the early universe. Interestingly, in a condensed matter 
context, the same correlation functions should enable us to estimate the superflow that 
would occur at a superfluid quench from fluctuations alonefTj. We shall make some steps 
in this direction. 

2 Defect Distributions 

Before going into details, some generalities about defect production will be useful. Our 
main interest is in vortices and our discussion will be centred about them. As will be 
seen, other defects are similar but simpler. The mechanism for vortex formation (termed 
the Kibble mechanism H in astroparticle physics) is well understood at a qualitative level. 
In this paper we shall only consider the simplest theory that permits vortices, that of a 
complex scalar field 0(x, t). The complex order parameter of the theory is {(f)) = rje 101 
and the theory possesses a global 0(2) symmetry that we take to be broken at its phase 
transition. Initially, we take the system to be in the symmetry-unbroken (disordered) 
phase, in which the field is distributed about = with zero mean. We assume that, at 
some time t = to, the 0(2) symmetry of the ground-state (vacuum) is broken by a rapid 
change in the environment inducing an explicit time-dependence in the field parameters. 
Once this quench is completed the 0-field potential takes the familiar symmetry-broken 
form V(4>) = -M 2 \<p 2 \ + A|0 4 |/4 with M 2 > 0. 

The transition for such a global symmetry is continuous and we expect that, as the 
complex scalar field begins to fall from the false ground-state into the true ground-state, 
different points on the ground-state manifold (the circle S 1 , labelled by the phase a of 
((f))) will be chosen at each point in space Q If this is so then continuity and single 
valuedness will sometimes force the field to remain in the false ground-state at = . 
For example, the phase of the field may change by an integer multiple of 2it on going 
round a loop in space. This requires at least one zero of the field within the loop, each of 
which has topological stability and characterises a vortex (or string). As to the density 
of the strings, if the phase a is correlated over a distance £, then the density of strings 

1 Had it been a first-order transition the field would have tunnelled towards the groundstate, with 
very different consequences for densities and correlations than those described here 
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passing through any surface will be 0(£~ 2 ) i.e. a fraction of a string per unit correlation 
area. The assumed lack of correlation of the field phase over larger distances than £ 
(inevitable on causal grounds in cosmological models) is interpreted as saying that the 
power P{k) of the field fluctuations, defined by 

(0(k)0*(k')) = (2Tt) 3 P(k)6 {3) (k - k') (2.1) 

has the form P(k) = 0(k°) for small k. That is, P(k) describes white noise at large 
distances. Several numerical simulations based on this assumption have been performed. 
The standard simulation, by Vachasparti and Vilenking], assumes a regular cubic lattice, 
in whose cells the 0(2) field phase is chosen at random. This assumption of 'white noise' 
leads to approximately 80% of the string network being in open string. While this 
number is lattice dependent || there is no doubt that a substantial fraction of string is 
open, or 'infinite'. We shall take these white noise field fluctuation predictions as a further 
benchmark against which to compare our dynamical predictions. On the completion of 
the transition, when fluctuations are too weak to eliminate or create strings, a network 
of strings survives whose further evolution is determined by classical considerations as 
the field gradients adjust to minimise the energy. 

Vortices are not the only defects permitted by global symmetries although, ultimately, 
they are the only ones of real interest to us. More generally, a global O(JV) scalar theory 
permits defects with integer topological charge in D = JV spatial dimensions (monopoles) 
or D = N + 1 spatial dimensions (vortices). Because of their relative complexity it is 
helpful to try out our methods on monopoles, which include the kinks on the line for a 
real scalar field theory (JV = D = 1) as a special case, and the monopoles in the plane 
for a complex scalar 0(2) theory (JV = D = 2). [The case of 0(3) global monopoles in 
three-dimensional space was examined by us elsewhere M and we shall not consider it 
further]. With minor qualifications^] a similar mechanism of continuous phase (or field) 
separation as that for vortex production is equally valid for the formation of these defects 
also. 

The question is, how can we infer these vortex, and other defect, densities and the 
density correlations from the microscopic field dynamics? The answer lies in the fact 
noted earlier, that for the global 0(2) theory of a complex scalar field 4> = (0i + <f>2) / V%, 
(0i, <p2 real) the string core is a line of zeroes of the fields <f) a (a— 1,2). The characterisation 
of an O(JV) global defect by its field zeroes is equally true for vortices in the plane (JV = 2) 
and kinks on the line (JV = 1). 

The problem is solved if we can identify those zeroes which will freeze out to define the 
late-time defects. This will require careful winnowing, since it is apparent that quantum 
fluctuations lead to zeroes of the fields on all distance scales (even in the disordered 
phase). For the moment we ignore this difficulty, and attempt to count every zero. The 
problem then reduces to that of determining the distribution of field zeroes, given the 

2 Because of the peculiarities of one spatial dimension, that would confuse the issue, we pretend that 
we are examining a one-dimensional section of a real field in higher dimensions and ignoring other degrees 
of freedom. 
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distribution of fields. This has been discussed in the literature on several occasions. We 
shall call repeatedly on the work of Halperin jlO| and Mazenko and Liu [11|, but see also 
Bray0. 



2.1 Kinks on the line 

As a prologue to the more difficult problem of vortices in three dimensions, we begin 
with the much simpler problem of identifying the zeroes of a real field 4>{x, t) in one 
f\ space dimension. To see how to proceed, consider an ensemble of systems evolving 
from one of a set of disordered states whose relative probabilities are known, to an 
ordered state as indicated above. At any given time t, the field will adopt one of the 
possible configurations 3>(sc), whose zeroes we wish to track. As the field evolves to its 
equilibrium values (one of the two minima of its potential V{<f>) = -M 2 2 + A0 4 / 4 ) th ese 
zeroes fluctuate and annihilate, but some of them will come to define the positions x of 
'kinks' (field interpolations from one minimum to the other at which <&'(x) > 0), some the 
position of 'antikinks' (at which $'(x) < 0), the one- dimensional counterparts of vortices 
and 'anti'-vortices. 

Suppose, at a given time, the zeroes of $(x) occur at x — Xi,x 2 , .... It is useful to 
define two densities. The first, 

p(x) ^Six-Xi), (2.2) 

i 

is the total density of zeroes, not distinguishing between kink zeroes and antikink zeroes 
(by which we now mean zeroes at which the field has positive or negative derivative). 
The second is the topological density, 

p(i)=^jiii(i-ii), (2.3) 

i 

where rij = sign($'(xj)), measuring (net) topological charge, the number of kink minus 
the number of antikink zeroes. 

Equivalently, in terms of the $-field, the total density is 

p(x)=5Mx)W(x)\, (2.4) 

since is the Jacobian of the transformation from zeroes to fields. Similarly, the 

topological density is 

p(x) = 6mx))&(x). (2.5) 

Analytically, it is not possible to keep track of individual transitions, but we can 
construct ensemble averages. If the phase change begins at time t then, for t > t , it is 
possible in principle to calculate the probability pt[&] that (j){x,t) takes the value 

3 See earlier footnote. 
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at time t. Ensemble averaging at time t is understood as averaging over the field 

probabilities Pt[$]. This is not thermal averaging since we are out of equilibrium. 

The situation we have in mind is one in which, for early times after the transition 
when the available space permits many domains, 

(p(x)) t = 0, (2.6) 

i. can equal likelihood of a kink zero or an antikink zero occurring in an infinitesimal 
length, compatible with an initially disordered state. However, the total zero density 

n(t) = (p(x)) t 

= fv§pt[§]5[$(x)]\&(x)\ >0 (2.7) 

is positive. The distribution of the zeroes is given by the density correlation function 

C(x;t) = (p(x)p(0)) t 

= Jv<$> pi[$]5[$(x)]o[$(0)]$'(x)$'(0), (2.8) 

{x 7^ 0) which will also be non-zero. 

To make this relationship more concrete, we observe that, on separating out the 



diagonal and non-diagonal terms in the expansion for p(x)p(y) from (|2.3|) , then 

p( x )p(y) = p(x)5(x - y) + g(x - y) (2.9) 

where 

g(x - y) = n i n j 5{x - Xi)5(y - Xj). (2.10) 

i¥=3 

That is, 

(p(x)p(0)) t = n(t)5(x) + C(x; t). (2.11) 
where C(x;t) = (g(x))t- Charge conservation 



dx(p(x)p(0)) t = (2.12) 

J — oo 

implies 

dxC(x;t) = -n(t), (2.13) 



requiring that C(x;t) = C(—x;t) be largely negative. 

We can now relate C(x; t) to the distribution and spacing of zeroes by calculating the 
variance of the topological charge 

n L = [ dxp(x) (2.14) 



o 
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on the interval I = [0, L\. In this particularly simple case til = — 1, 0, 1. Then 

(A t n L ) 2 = (nl) t 

rL rL 

= dx dy (p{x)p(y)) t 
Jo Jo 

= Lh(t) + / dx dyC(x — y;t) 
Jo Jo 



L dx [ L dyC(x-y;t). (2.15) 



x<0 



from (ET 



If h(x;t) is defined by C(x;t) = dh(x;t)/dx then 

d(A t n L f 



dL 



-2h(L;t). (2.16) 



However, if p^ is the probability that, on average, a length L of the line contains an odd 
number of zeroes, then (A t n^) 2 = p^. Thus 

C(L;t) = h'(L;t) = - 1 -^. (2.17) 

As an extreme case we note that, for an array of equally spaced zeroes, separation £(£), 
Pl is saw-toothed, period 2£, from which 

h(x; t) = ^^ysign(sin(7rx/e(t))), (2.18) 

whence C(x;t) = h'(x;t) is a sum of ^-functions. At the other extreme, independent 
(Poisson) zeroes, mean separation £(t), give 

h(x;t) = e-W*® /2£(t) (2.19) 

and 

C{L-t) = -e- 2L/m /i 2 {t). (2.20) 

These provide a useful guide, and we shall return to them later. However, some 
caution is necessary since, on the line, kink zero is followed by anti-kink zero and vice- 
versa. As an intermediate step between this and the vortices of the complex field 0(x, t) in 
three dimensions, we shall consider the same field in two dimensions, for which monopoles 
(point defects) occur. 



2.2 Monopoles in two dimensions 

An 0(2) complex field 3>(x) in the plane permits monopoles, at this stage identified by 
zeroes of the field with non-trivial winding number. After the transition is completed, 
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the relevant zeroes will define the centres of effectively classical monopoles. Although 
Derrick's theorem prohibits finite-energy static monopole solutions to the classical field 
equations, a non-zero monopole density provides a cut-off to the logarithmic tails of the 
individual monopoles, and there is no problem on this score. 

For the moment we continue to count every field zero, at whatever scale and however 
transient. If <&(x) has zeroes at X!,x 2 ,X3, ... the total and topological densities p(x) and 
p(x) are the straightforward generalisations of ( |2.2|) and ( [2. 3D , 

p(x) = 5>(x- Xi ), (2.21) 

i 

and 

p(x) =5>^(x- Xi ), (2.22) 

i 

where rti = ±1 is the winding number of the zero (higher winding numbers are taken as 
multiple zeroes). 

The relationship between zeroes and fields is through the Jacobian, giving 

p(x) = 5 2 [$(x)] ejfc ^$ 1 (x)^<l> 2 (x), ij = 1,2 (2.23) 

where eu = —£21 = 1, otherwise zero. As before we assume that it is possible in principle 
to calculate the probability pt[$] that 0(x, t) takes the value $(x) at time t from which 
we define the ensemble averages (_F[$]) t . The topological charge average (p(x)) t is taken 
to be zero as before. The non-zero total density and the topological density correlation 
functions are defined as for the one-dimensional case, 

n{t) = (p(x)) t 

= J V$ pti^M^kd^Wdk^Wl >0 (2.24) 

and 

C(r;f) = (p(x)p(0)) 4 

= Jv® p t [$]5 2 [$(x)]5 2 [$^ 

(r = |x| 7^ 0). Charge conservation again applies, as 

J d 2 xC(r;t) = -n(t), (2.26) 

but it is not so easy to give a direct interpretation to C(x;t) in terms of monopole 
separations. 

4 Throughout, it will be convenient to decompose $ into real and imaginary parts as $ = -^( < J>i+« ( I?2)- 
This is because we wish to track the field as it falls from the unstable ground-state hump at the centre 
of the potential to the ground-state manifold in Cartesian field space. 
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Nonetheless, it is still useful to consider the variance in the topological charge 

n s = I rf 2 xp(x), (2.27) 
Jxes 

this time through a region S in the plane, (area s) since it can be an observable quantity. 
For example, on quenching 4 He in an annulus S, (A t ns) 2 is a measure of the supercurrent 
generated by the quench [0] . Nothing that we have said so far really requires a relativistic 
theory, so let us pursue it a little further. Then, from ( |2.26|) it follows that 

(A^ 5 ) 2 = / d 2 xf d 2 y(p( X )p(y)) t 
Jxes Jyes 

= - / d 2 x I d 2 yC{\x-y\;t) (2.28) 

Jx^S JyeS 

the two-dimensional counterpart to Q2.15|) . Suppose that there are short-range vortex- 
zero antivortex-zero correlations in which C(r;t) is short-range in r, with lengthscale 
f(t). If f(t) is the only lengthscale then C(r;t) is 0(£~ 4 (t)). 

With x outside S, and y inside S, all the contribution to (A t n,s) 2 comes from the 
vicinity of the boundary of S, rather than the whole area. More precisely, suppose that 
S is a disc of radius L £]. Then ( |2.28j ) can be written as 

pL /"oo flit 

(A t n s ) 2 = -2n rdr r' dr' dO C(R;t) (2.29) 
Jo Jl Jo 

where 

R 2 = r 2 + r' 2 - 2rr'cos6. (2.30) 

The integration region in r,r' is restricted to \r' — L\ = 0(£(t)) — \L — r\. If C(r;i) is 
non-singular at r = and not varying too rapidly at the origin, then the ^-integration is 
limited to a range 0(£(t)/L), whence 

(A t n s ) 2 = 0(^J. (2.31) 

This robust result is compatible with, but does not imply, a random walk in field phase 
around the perimeter and we shall return to it later, when we shall see that our model 
implies short-range correlations 0. 

2.3 Vortices in Three Dimensions 

Finally reaching our goal of vortices in three dimensions, we define the topological line 
density of zeroes p(r) [[K], |ll|] by 

r dTl 

p(x) = £ J ds^5 3 [x - R n (s)]. (2.32) 

J The extension to an annulus is straightforward. 

6 At the other extreme we note that, if the monopoles and antimonopoles are individually and mutually 
uncorrelated (a double Poisson distribution), then (Ans) 2 = sn. That is, the variance grows linearly 
with area. 



In (2.1) ds is the incremental length along the line of zeroes R n (s) (n—1,2,.. .) and ^ L 
is a unit vector pointing in the direction which corresponds to positive winding number. 
As in the previous case, we begin by counting all zeroes, anticipating that the relevant 
ones will become the cores of the strings of the resulting network. Yet again, a finite 
string density once the transition is complete provides the cutoff necessary to eliminate 
the infrared divergent tails of the energy densities that we would expect from Derrick's 
theorem. 

It follows that, in terms of the zeroes of $(x), p»(x) can be written as 

Pi (x) = 5 2 [$(x)] e ^^$ 1 (x)9 fc $ 2 (x), (2.33) 

where <5 2 [<I>(x)] = 5[$i(x)]<5[$ 2 (x)]. The coefficient of the 5-function in ( p73| ) is the Jaco- 
bian of the more complicated transformation from line zeroes to field zeroes. We shall 
also need the total line density p(x), the counterpart of p(x) of 



Pi(x) = 5 2 [$(x)]|e^^$ 1 (x)9 fc $ 2 (x)|. (2.34) 

As before, ensemble averaging (F[§]) t at time t means averaging over the field prob- 
abilities . Again we assume 

(Pi(x))* = 0. (2.35) 

i. e. an equal likelihood of a string line- zero or an antistring line-zero passing through an 
infinitesimal area. However, 

n(t) = (p i (x)) i > (2.36) 

and measures the total line-zero density in the direction i, without regard to string 
orientation. The isotropy of the initial state guarantees that n(t) is independent of the 
direction i. Further, the line density correlation functions 

C ii (x;*) = <Pi(x)p i (0)) t (2.37) 

will be non-zero, and give information on the persistence length of line zeroes. It will be 
convenient, for later work, to decompose Cy(x;£) as (r = |x|), 

1 3 \ id ( -l\ { " Xj t^ J 3 



C i3 (x; t) = A(r; t) ( 5 i3 - \ + B(r; t) (=^2 ) , (2.38) 



The realisations of n(t) and CV,(x; t) in terms of pt[<f>] are simple generalisations of ( 2.24 ) 
and ( |2.25| ) and will not be given explicitly. 
Charge conservation now means 

J d 3 x C i3 -(x;t) = (2.39) 

without any inhomogeneous term. On taking the trace in Q2.39p it follows that 

d 3 x (2A(r; t) + B(r; t)^j = (2.40) 
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We note that, from (|2.32 ), the integral of Pj(x) over an open surface S with normal in 



the j-direction does not measure the winding number along the boundary of S since each 
line zero is weighted by the cosine of its angle of incidence on S. Thus the variance of 
the winding number (topological charge) through S is essentially the two-dimensional 
problem discussed previously. 

There are no simple quantities that can be calculated for line zeroes (e.g.Poisson 
strings). However, we note that, if x = (0,0, r), then 

di(x;t) = C 22 (*;t) = A(r;t) (2.41) 

measures the likelihood of there being a string zero (or antistring zero) in the same 
direction at separation r. A negative value of A of 0(n 2 ) indicates the presence of 
vortex- antivortex line zero pairs at separation r. Further, for the same r, 

C 33 (x;t) = B(r;t) (2.42) 

measures the tendency of a vortex to bend on a distance r (provided r is sufficiently 
small for it to be measuring the same string). Once Bn 2 is negligable the line has bent 
away from its initial direction. 



3 Ensemble Averaging 

In relating distributions of defects to field fluctuations we have seen that we need the 
field probability at all times. There is no difficulty in writing a formal expression 
for pt[$], although its calculation is another matter. Details are given in [|J, and we shall 
only provide the briefest recapitulation. 

Take t = t as our starting time. Suppose that, at to, the system is in a pure state, 
in which the measurement of (j) would give $o( x ). That is:- 

0(to,x)|$ o ,to) = $ |$o,*o>- ( 3 - 43 ) 

The probability pt . [<&/] that, at time tf > to, the measurement of will give the value $/ 
is p tf [§f] = \^fo\ 2 , where \I//o is the state- functional with the specified initial condition. 
As a path integral 

*/o= / V<j>exv\iS t [<j>}\, (3.44) 

where S t [(j)] is the (time-dependent) action that desribes how the field <fi is driven by the 
environment, V<p = Yla=i an d spatial labels have been suppressed. Specifically, for 
t > to the action for the field is taken to be 

St[<f>] = J d D+1 x(±d^ a d^ a - \m 2 (t)<p 2 a - \\(tM) 2 y (3.45) 

where m(t), X(t) describe the evolution of the parameters of the theory under external 
influences, to which the field responds. The spatial dimension D = 3 for the relevant 
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case of vortices, and D = 2 for monopoles in the plane []. As with it is convenient to 
decompose the complex field i n terms of two massive real scalar fields a , a — 1, 2 as 
= (0i + i<t>2) I v2) i n terms of which S[(f>] shows a global 0(2) invariance, broken by the 
mass term if m 2 (t) is negative. 

It follows that pt f [$f] can be written in the closed time-path form 

Pt f [$ f ] = / 0±(t/)_$/ exp{^[0 + ] - S t y>.])}. (3.46) 

Instead of separately integrating 0± along the time paths to <t <tf, the integral can be 
interpreted as time-ordering of a field along the closed path C + © C_ where = + on 
C + and = 0_ on C_. The two-field notation is misleading in that it suggests that the 
+ and 0_ fields are decoupled. That this is not so follows immediately from the fact 
that 0+(t/) = 0-(t/). It is necessary to keep this in mind when we extend the contour 
from tf to t = oo. Either + or 0_ is an equally good candidate for the physical field, 
but we choose + . With this choice and suitable normalisation, p t becomes 



Figure 1: The closed timepath contour C + © C_. 



Pt f [$f] = I D0 + X?0-5[0+(t) - $/] exp{z(^ t [0 + ] - $ [<£_])}, (3.47) 

where 5[4> + {t) — $/] is a delta functional, imposing the constraint + (t, x) = $/(x) for 
each x. 

The choice of a pure state at time to is too simple to be of any use. The one fixed 
condition is that we begin in a symmetric state with (0) = at time t = to- Otherwise, 
our ignorance is parametrised in the probability distribution that at time t , 0(t o ,x) = 
$(x). If we allow for an initial probability distribution Pt [$] then p tf [§f] is generalised 
to 

Pt f [$ f ]= j-DQptM J V<f> + V<f>_6[<t> + {tf)-®f} exp{z^[0 + ]-^[0_])}. (3.48) 

7 D = 1 for kinks on the line 
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It is impossible to derive p t analytically for general initial conditions. Fortunately we 
shall see that, in many circumstances, the details of the initial condition are largely irrel- 
evant. All the cases that we might wish to consider are encompassed in the assumption 
that $ is Boltzmann distributed at time to at an effective temperature of To = (3q 1 ac- 
cording to a Hamiltonian i?o[^*], where the subscript zero denotes to rather than a free 
field. That is 

p t M = ($,to|e-^|$,to) = / V<f> 3 exp{iS [<f> 3 }\, (3.49) 

for a corresponding action So[03], in which 3 is taken to be periodic in imaginary time 
with period (3 . We take S [(j) 3 ] to have the standard form in 3 as 



d D+1 



X 



(3.50) 



We stress that m , Ao and (3q parametrise our uncertainty in the initial conditions. The 
choice (3q — > oo corresponds to choosing the ptW\ to be determined by the ground state 
functional of H , for example. For the sake of argument we take T = Pq 1 to be a 
temperature higher than the transition temperature T c and m = m(T ) {ni^ > 0) to be 
the effective mass at this at this temperature. Whatever, the effect is to give an action 
Ss[cf>] in which we are in thermal equilibrium for t < t during which period the mass 
m(t) and coupling constant A(t) take the constant values m and Ao respectively. 
We now have the explicit form for p tf [$ /] , 

p tf [$ f ] = J B V<f> 3 V<f> + V<f>_ exp|^ o [03] +i(S[<j> + ] - 5[0_])}5[0 + (t 7 ) - (3.51) 

where the boundary condition B is 0±(to) = 03 (to) — ^3(^0 — iPo)- This can be written 
as the time ordering of a single field:- 

= I Vct>e iS ^ 5[Mtf) ~ <*>/], (3.52) 

along the contour C = C + © C_ © C 3 , extended to include a third imaginary leg, where 
4> takes the values <f> + , 0_ and <p 3 on C+, C_ and C 3 respectively, for which Sc is 5'[0+], 
S*[0_] and (So [03]. Henceforth we drop the suffix / on $j and take the origin in time from 
which the evolution begins as t = 0. 

We perform one final manoeuvre with p t [$] before resorting to further approximation 
to demonstrate how we can average without having to calculate explicitly. Further, 
this will enable us to avoid a nominally ill-defined inversion of a two-point function later 
on, a consequence of the seeming independence of + and 0_ mentioned earlier. Consider 
the generating functional:- 

Z\j+,j-,h] = J B V(j> exp|^ c [0] +i J j'0|, (3.53) 



12 



to 



C 3 

to - ifio 



Q ; m(r) 



c 



5Re(r) 



Figure 2: A third imaginary leg 
where / j(p is a short notation for:- 

/foc 
3<j> = J Q dt [j + (t)</>+(t) - j-<j>-(t)] + J o Mt)Mt)dt, (3.54) 

omitting spatial arguments. Then introducing a (x) where a = 1, . . . , N, we find:- 

p t/ [$] = JvaJ V<j> exp|^ c [0]J expj* | rf D xa a (x)[0 + (t /? x) - $(x)] a J 

= Jva expj-i ^ a a $ a j Z[a, 0, 0], 

where a is the source a(t, x) = a(x)5(t—tf). As with £>0, Pa denotes Y\i T^o: a . Ensemble 
averages are now expressible in terms of Z^. Of particular relevance, W a b(|x — x'\;t) = 
($ a (x)$ 6 (x / )) t is given by 

W^(|x-x'|;t) = Jv$ $ a (x)$ 6 (x') J Va expj-* J a a $ a | Z^a, 0, 0] 



Va 

On integrating by parts 

W ab (\x-x'\;t) 



5 2 



5a a (x.)5ab(x.- 
(0a(x,t)0 fe (x',t)), 



;r^[a,0,0] 



a=0 



(3.56) 



the equal-time thermal Wightman function with the given thermal boundary conditions. 
Because of the non-equilibrium time evolution there is no time translation invariance in 
the double time label. 
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4 A Gaussian Model for Defect Distributions 



We have yet to specify the nature of the quench but it is already apparent that, if we 
are to make further progress, additional approximations are necessary. We return to our 
one-dimensional example. 

4.1 Kinks 

Suppose the fields and their derivatives at the same point show approximate indepen- 
dence. Then the zero density n(t) of Q2.7Q separates as 

fiW^i^)])^)!),. (4.57) 

and we can estimate, or bound, each factor separately. Specifically, on writing the 5- 
function term as 

(6[$(x)]) t = [ ^a(exp{i / dyj(y)<S>(y)}) t (4.58) 



(4.59) 



where j(y) = a5(y — x), it follows that 



n{t) < 



1 


W"(0; t) 


V2vr 


W(0;t) 



1 + Connected Correlation Functions 



on using the Schwarz inequality, where primes on W denote differentiation with respect 
to x. Some care is needed. The connected correlation functions that appear in ( 4.59) 



are 0(A) but depend on time, growing as the field grows away from = 0. As is well 
understood, perturbation theory breaks down at long enough times. However, for small 
times and weak coupling the first term in (|4.59|) is reliable. We see that, if the Fourier 
transform W(k;t) of W(x;t) is dominated by wave vectors k(t) = 0(£ _1 (£)) at time t, 
then 

n(t)<o(^). (4.60) 

This shows that sets the domain size. 

A similar qualitative analysis can be attempted for the density correlation functions, 
but with less success, in the absence of further approximation. To be concrete, let us 
consider the consequences of pt[3>] being Gaussian, for which the approximate equality 
in ( 4.57|) becomes exact. That this is not a frivolous exercise in solving what we can 



solve but a representation of reasonable dynamics will be shown, in part, later, for the 
'slow-roll' dynamics that we shall adopt. For the moment we take it for granted. 

Specifically, suppose that (still for the one- dimensional case) $ is a Gaussian field for 
which 

($(x)) t = = ($(x)&(x)) t , (4.61) 

and that 

(<S>(x)$(y)) t = W(\x-y\;t). (4.62) 
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All other connected correlation functions are taken to be zero. Then all ensemble averages 
are given in terms of W(r; t) which we have have seen to be the equal-time Wightman 
function, 

W{\x-y\;t) = {<f>{x,t)<t>{y,t)) (4.63) 

with the given initial conditions. In our case, where we shall assume thermal equilibrium 
initially, this is the usual thermal Wightman function. It is then straightforward to see 
that, if 

then 

n(t) = -(-f(0;t))l (4.65) 

7T 

in agreement with (|4.59|) . On using the same exponentiation that (<5[$(x)])i equals 
/ tf,a{e ia ^ x ') t it takes only a little manipulation to cast the correlation function C(\x\;t) 
of ( j2.8|) in the form 

G(r;() = ^M, (4.66) 

where 

Kr;t)= ~ f(r;t) ■ (4.67) 
2ny/l-f*(r;t) 

If we are given C(r;t), and wish to infer f(r;t), and hence its power spectrum, we 
integrate equation ( |4.67| ) as 



/(r; t) = sin 



- - 2tt / dr'h(r';t) 



(4.66 



in which h(r;t) is, in turn, derived from C(r;t) through equation (|4.66 ). For example, 
for the Poisson distribution of zeroes given earlier in fl2.19|) , we find 

f(r;t) = stn(^e- 2r/m ^j (4.69) 

However, the situation is usually the converse, with the dynamical model predicting 
f(r;t), from which C(r;t) follows. As a concrete example, suppose that 

W{r;t) = J 4ke ikx W{k;t) (4.70) 

where 

W{k-t) = 5{k 2 -k&t)). (4.71) 
That is, there is only one wavelength in the model. Then 

f(r;t) = cos (k (t)r), 

n(t) = k (t)/n. (4.72) 
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The resulting 

h(r;t) = ^-sign(sm (k (t)r)), (4.73) 
is exactly that of ( [2.18 ), in which the correlation function C(r;t) is a sum of delta- 



functions, corresponding to a regular array of zeroes. 

Now suppose that, instead of a single wavelength, W(r; t) has contributions from 
frequencies peaked about k Q , as 

W(k;t) = I^xp{-i(A;-A; (t)) 2 /(Afc (t)) 2 } + exp{-i(A; + fco(t)) 2 /(AA;o(t)) 2 }) (4-74) 
In the limit Ak Q (t) —>■ we recover the regular array but, for Ako(t) ^ 0, 

W{r-t) = J dk coskx exp{-i(A;- k (t)) 2 / (Ak (t)) 2 } (4.75) 

giving 

f(r; t) = (cos(k (t)r)) exp{-K 2 (Ak (t)) 2 }. (4.76) 
The corresponding h(r; t) is no longer that of ( |2.18| ), but 

h(r) = feo^(V)+r(Afc ) 2 cos(fc r) exp{ _ l r2(Afco)2} , (4J?) 
— (cos 2 (k r)) exp{—r 2 (Ak ) 2 } 2 



and we have dropped the time-labelling. A cursory inspection shows that h(r) is damped 
periodic, and the broadening of its derivative C(r; t) from its 5-function behaviour when 
Afco(t) 7^ measures the variance A£(t) in the separation £(£) between adjacent zeroes. 
We see that, for small Ako(t) / ko(t) , 

with a coefficient of proportionality of 0(l)f]. We note that, even when the variance 
A£(t)/£(t) is large, so that the zeroes appear much more random, the Gaussian behaviour 
is very different from the linear exponential behaviour of the Poisson distribution noted 
earlier. There is always correlation. 

4.2 Monopoles 

The extension to 0(2) line zeroes is messy, but leads to no surprises. Specifically, suppose 
that 

($«(x)) t = = (^ a (x)d^ b {x)) t , (4.79) 
8 It follows that the effect of the variance is to increase the density of zeroes to n{t) = 

y/k*(t) + (Afc (i)) 2 A- 
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and, further, that 

($ a (x)$ 6 (x')) t = W ab (\x - x'|; t) = 5 ab W(\x - x'|; t), (4.80) 
is diagonal. As before, all other connected correlation functions are taken to be zero. 



The density calculation proceeds as before. It follows [TO, [TT[ that 



n{t) 



2tt 



-f"(0;t)). 



(4.81) 



where derivatives are taken with respect to r = |x — x'|. The zero density- density 
correlation function is more complicated than for the kinks, albeit still in terms of h(r; t) 

of qnj7p , as|n 

C(r;t) = -h(r;t)h'{r;t). 



If C(r; t) is given, we can integrate equation (|4.82|) to obtain f(r; t) as 



f(r;t) = sin 



7T 



2tt / dr 



dr' r' C(r'; t) 



1/2 



(4.82) 



(4.83) 



on incorporating the conservation of topological charge. 

However, we are more interested in calculating C(r; t) for specific f(r; t). As for kinks 
above, the simplest case corresponds to taking 

(4.84) 

(4.85) 

(4.86) 
(4.87) 

(4.88) 



W(r;t) = J d 2 ke lk ~W(k;t) 

in which there is only one wavelength 

W(k ] t) = 6(k 2 -k 2 (t))- 
where k (t) depends on time t. As a result 

f{r;t) = J (k (t)r), 



n{t) 



71 



and, from (|4.67|) , 



h(r;t) 



Ji{k {t)r) 



2tt^1 - J 2 (A; (t)r) 

Because of its rotational invariance C(r; t) of ( 2.25|) cannot show the 5-function behaviour 
of its counterpart on the line, but regularity is implicit in the strong regular oscillatory 
peaking of C that comes from the Bessel functions in the numerator. 

If, instead of a single wavelength, W(r; t) has contributions from frequencies peaked 
about ko with variance Ako, then the secondary peaking is reduced. Specifically, consider 



W(r;t) 



dk J (kr) exp{--(k - k (t)) 2 /(Ak (t)) 2 }. 



(4.89) 
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In fact, unlike the one-dimensional case, for which W(r; t) and f(r; t) of ( |4.76| ) retain their 
oscillatory behaviour, when Ak (t) / k (t) is large enough W(r; t) does not oscillate or even 
change sign. This is understood as indicating a less regular distribution of monopoles. 
The dispersion in wavelength Ako(t) leads to a variance A£ in their separation that we 
expect to satisfy 

M^^M (490) 

In itself this short-range behaviour of f(r; t) is enough to give the variance in winding 
number through a surface of perimeter L as (AiV) 2 = O(L) for the reasons given earlier. 



4.3 Vortices 

Finally, for line zeroes the Gaussian approximation has yet more complicated conse- 
quences. As for monopoles, we assume diagonal correlation functions 



($ a (x)$ 6 (x')) t = W ab (|x-x'|;t) = 6 ab W(\x-x.'\;t), 



(4.91) 



from which /(r; t) can be defined as in (|4.64 ), and vanishing field expectation value and 
and the independence of the field and its derivatives 



($a(x))t = = ($ a (x)^$ 6 (x))i, 

As might have been anticipated, n(t) is as for monopoles 



(4.92) 



n(t) = -(-f'(0;t)), 

Z7T 



(4.93) 



whereas flOl, PH the transverse and longitudinal parts of the density correlation ( |2.38|) 
are, still in terms of h(r; t), 

2 

A(r;t) = -h(r;t)h'(r;t) (4.94) 
(the same in form as C(r; t) of (|2.25|) ), and 



B(r;t) 



T /i 2 (r;t) > 0. 



(4.95) 



We note that B is positive. The conservation law (|2.40|) follows, as 

d 



J d 3 x (2A 



B) oc / dr — ( rh 2 
lo dr 







(4.96) 



Given either A(r; t) or B(r; t) we can reconstruct /(r; t). For given A(r; t), f(r; t) has 
a similar form to that of (|4. 83|) (except for the absence of charge conservation). More 
simply, for given B(r;t), 



f(r;t) = sin 



--nV2 I" dr'r'(2B(r';t)) 1 / 2 
2 Jo 



(4.97) 
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If, as before, we asssume a single wavelength, then 

W(r;t) = J fke^5(k 2 -k 2 (t)) 
sin(/c Q (t)r) 



ko(t)r 



sine (fco(t)r) (4.98) 



It follows that 



and 



nit) = ^ (4.99) 
w ^ = sin(fc (t)r) - fc (t)r cos(fc (t)r) ^ 



2nr^k 2 (t)r 2 - sm 2 (k (t)r) 

from which A(r; t) and f?(r; £) can be constructed from above. A and B again show peri- 
odic peaking, pointing to as regular a distribution of line-zeroes as possible. Specifically, 
when k (t)r ^> 1 

A[r . t) ^zMl^p^ (4101) 



and 

cos 2 k (t)r 
2irr 



up to non-leading terms. On the other hand, when ko(t)r <C 1 A(r;t) = 0{n 2 {t)) and 
negative, whereas B(r;t) = 0{n 2 (t)/k 2 {t)r 2 ). 

Yet again, if, instead of a single wavelength, W(r;t) has contributions from frequencies 
peaked about k with variance Ak (t) the secondary peaking is reduced. If we take 

W(r;t) = Jdksinc(kr) exp{-^(A;- k (t)) 2 / (Ak (t)) 2 } . (4.103) 

then, as for monopole zeroes in two dimensions, once Ako(t)/ko(t) is sufficiently large, 
W(r;t) ceases to oscillate or even to vanish. Yet again we understand the variance in k 
as inducing a variance A£ in their separation approximately satisfying ( |4.9U| ). However, 
the short- distance behaviour of A and B is unchanged qualitatively. 

We shall see later, in our model making, how distributions approximately of the form 
( [4.103] ) arise naturally. However, there is another, even simpler, and not wholly dissimilar 



way of introducing a dispersion about the wavenumber ko(t) that has been adopted in 
a cosmological context 0. Because of the similarities between it and the model that we 
shall introduce later we shall look at it in some detail. The work of essentially consists 
in taking W(r; t) as 

, , rko(t) / k \ 2+n 

W (n) (r;t)oc dksmc(kr) [-—) (4.104) 

Jo \k (t)J 

describing fields with a variable power spectrum k n , cut off at k = k Q (t). 
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For n > the behaviour of the normalised field correlation functions f( n \r;t) = 
W^ n \r\ t)/W^ n \0; t) is determined by the sharp cutoff, unlikely to be present in any 
realistic models. However, it follows that the limit n — > oo reproduces the single-mode 
result of ( |4.98| ). The density n(t) is calculated easily as 

reducing to ( |4.99| ) as n — > oo, as it must. 

As n decreases the power in long wavelength modes increases. For example, for 
n = 0,-1,-2 

/ (0) M) = 77. 7I^ ( sin ( fe o(t)r) - ko(t)r cos(fc (t)r) 



(ko(t)r] 



/ M) M) = --^(l - cos(k (t)r)) (4.106) 

f ( - 2) (r;t) = -J— Si{ko(t)r) 
ka(t)r 

The case n = corresponds to white noise on scales larger than fc^ 1 ^. We see immediately 
that f^~ l \r] t) of fl4.106| ) has double the period of the large-n W(r; t) of ( |4.98| ), indicating 
a reduced density directly. 

If we calculate the density-density anticorrelation function A(r; t) for these and 
others we find that, for small r (where the effects of the sharp cutoff are less apparent) 
it can be written as 

A(r;t) « -n 2 (t)(a n - b n (k (t)r) 2 ) (4.107) 

where a n and b n are positive and, by extracting a factor n 2 (t) we are working at constant 
density. The strength a n of the anticorrelation is approximately constant and 0(1) for 
n > 0, whereas for n < it diminishes rapidly, vanishing when n rs —2. For n > — 2 
the range of the anticorrelation is determined by b~ x . Over the range —2 < n < oo, b„ 
is approximately linear, growing from approximately zero at n = —2 with slope 0(1). 
Thus, as n becomes increasingly negative (but n > —2) the range over which the strings 
have influence on one another increases, but with diminishing strength. For n = —2, 
ko(t) ceases to be a dominant wavenumber and there is no quantity to identify with a 
domain size. We would interpret this as implying the greatest variance in string-zero 
separations. 

For large r, we note that 

/( " ,(ri ') = °((w) (4108) 

for integer n > —1, and f(~ 2 \r;t) = O(l/(ko(t)r)). The vanishing of n(t) at n = —3 
is a consequence of the infrared divergence at this value, but can be ignored since it is 
difficult to see how values of n < — 2 could arise from realistic dynamics. 

9 White noise on all scales gives a (5-function for W. 
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With these examples behind us we are now ready to attempt to make predictions in 
a simple dynamical model of Gaussian fluctuations. 

5 Gaussian Dynamics from an Instantaneous Quench 
and its Coarsening 

In practice, our ability to construct pt[<&] or, equivalently, the field correlation functions 
over the whole timescale t > to from initial quantum fluctuations to late time classicality 
is severely limited. In particular, a Gaussian approximation can have only a limited 
applicability. To see what this is, it is convenient to divide time into four intervals, to 
each of which we adopt a different approach. With M setting the mass-scale, there is 
an initial period to < t < U — (^(M^ 1 ), before which the field is able to respond to the 
quench, however rapidly it is implemented, and which we can largely ignore. Assuming 
weak coupling, of which more later, this is followed by an interval U < t < t sp in which, 
provided the quench is sufficiently rapid, domains in field phase form and grow. For 
the reasons given above vortices (or monopoles) will appear and be driven apart by the 
coalescence of these domains. For the quartic winebottle potential V((f>) with minima 
at |0| =0o domain growth begins to stop once the field has reached its spinodal value 
|0| 2 = 0q/3 = 0(M 2 /A) at the ring of inflection V"((p) = 0. This occurs at times peaked 
around some t = t sp . To estimate t sp we observe that, for short times and weak coupling, 
the long wavelength modes with wavenumber k m grow as 0k(O — 0(Me Mt ) as they 
fall from the top of the potential hill. Thus t sp is given in terms of the coupling strength 
by 

exp{-2Mt sp } = 0(A) (5.109) 

or equivalently, t sp = 0(M~ l ln(l/A)). We assume that A is sufficiently small that t sp is 
significantly larger than ij. The defects are beginning to freeze in and fluctuations are 
now too small to undo them. In the third period, beginning at t sp , the field magnitude 
relaxes dissipatively to the ground state values and the vortices complete their freezing 
in. Finally, in the last period, the vortices behave semiclassically. 

We shall see that our ability to calculate from first principles is limited to the second 
period tt < t < t sp . On the (as yet unproven) assumption that the distribution of relevant 
zeroes at time t sp is left largely unchanged by their final freezing in, this distribution of 
vortices can then be taken as initial data for the final evolution of the network. Similar 
considerations apply to monopoles. 

It is not difficult to justify our adoption of Gaussian field fluctuations for this second 
period U < t < t sp of vortex production. We have already assumed that the initial 
conditions correspond to a disordered state. In the absence of any compelling evidence 
to the contrary we achieve this by adopting thermal equilibrium at a temperature T 
higher than the critical temperature T c for t < t , as we anticipated earlier. That is, the 
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action So of ( |3.50| ) that characterises the initial conditions is 



Sol 



d 



D+l 



X 



(5.110) 



with m 2 (T) > 0. If Ao is weak the resulting field distribution is approximately Gaussian, 
and it is little loss to take it to be exactly Gaussian, Ao = 0. As we shall see later, initial 
conditions generally give slowly varying behaviour in the correlation function, in contrast 
to the rapid variation due to the subsequent dynamics, and calculations are insensitive 
to them. 

In order to have as simple a change of phase as possible, we assume an idealised 
quench, in which, at t — to, m 2 (t) changes sign everywhere. Most simply, this change in 
sign in m 2 (t) can be interpreted as due to a reduction in temperature. Even more, at first 
we further simplify our calculation by assuming that, for t > t , m 2 (t) takes the negative 
value m 2 (t) = —M 2 < immediately, where — M 2 is the mass parameter of the (cold) 
relativistic Lagrangian. That is, the potential at the origin has been instantaneously 
inverted, breaking the global 0(2) symmetry. If, as we shall further assume, the A|0| 4 
field coupling is very weak then, for times M(t — t ) < ln(l/A), the 0-field, falling down 
the hill away from the metastable vacuum, will not yet have experienced the upturn of 
the potential, before the point of inflection and we can set A = 0. Thus, for these small 
time intervals, pt[&] is Gaussian, as required. Henceforth, we take to = 0. 

For such a weakly coupled field the onset of the phase transition at time t = is 
characterised by the instabilities of long wavelength fluctuations permitting the growth 
of correlations. Although the initial value of ((f)) over any volume is zero, we anticipate 
that the resulting evolution will lead to domains of constant <j) phase, whose boundaries 
will trap vortices. This situation of inverted harmonic oscillators was studied many years 
ago by Guth and Pi|l3|] and Weinberg and Wu||14j|. In the context of domain formation, 
we refer to the recent work of Boyanovsky et al.|16|], and our own[|J. 

Consider small amplitude fluctuations of <p a , at the top of the parabolic potential 



hill. Long wavelength fluctuations, for which 
If their growth rate Q(k) = 



|k| < M , begin to grow exponentially. 
I M 2 — |k| 2 is much slower than the rate of change of 
the environment, then those long wavelength modes are unable to track the quench. 
Unsurprisingly, the time-scale at which domains appear in this instantaneous quench is 
ti = 0(M~ 1 ). As long as the time taken to implement the quench is comparable to ti 
and less than t sp = 0(M~ l ln(l/A)) the approximation is relevant. We shall relax the 
condition in the next section. 

For the moment we ignore the effect of the interactions that stabilise the potential 
and permit the zeroes corresponding to the defects to freeze in. In this free-roll period, 
W(r;t) has to be built from the modes U^ k , 
motion 

" d 2 



U^ k = {U ak )*, satisfying the equations of 



,11- 



+ k 2 + m 2 (t) 



(5.111) 
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For the idealised case proposed above m 2 (t) is of the form 



m 2 {t) = m 2 > if t < 0, 
= —M 2 < if t > 0. 



(5.112) 



If we make a separation into the unstable long wavelength modes, for which |k| < M, 
and the short wavelength modes |k| > M, then W(r;t) is the real quantity 



W(r;t) 



k^c{k)uut)u- k {t) 



\V\<M 



4 u ke lk -*C(k) 



+ 



k^C{k) 



'|k|>Af 

with no summation over a, r = |x| and 

n 2 (k) = 

w 2 (k) = 
A(k) = 

a{k) = 



1 + A(k)(cos\i(2VL(k)t) - 1 
l + a(Jfe)(cos(2u;(Jfe)*) - 1 



(5.113) 



(5.114) 



iir - |k| 2 , 

—M 2 + |k| 

1 / UJ' 

2\ 1 + n 2 (k)J' 

2 V w 2 (k) 



(5.115) 



The initial conditions are encoded in C(k), which takes the familiar form 

1 



C(k) 



2u(k) 



coth(l3 uj(k)/2) 



(5.116) 



in which uj 2 (k) = |k| 2 + m\. 

In the single-mode and other approximations adopted earlier the folly of counting 
all zeroes was not apparent. It is now, in the presence of the ultraviolet divergence of 
W (r; t) at r — in all dimensions. None of the expressions given above is well-defined. 
To identify which zeroes will turn into our vortex network requires coarse- graining. 

The way to do this is determined by the dynamics. Firstly, we note that if $ is 
Gaussian, then so is the coarsegrained field on scale L, 



$ L (x) = J dV/Gx-x'Q^x') 



(5.117) 



where J(r) is an indicator (window) function, normalised to unity, which falls off rapidly 
for r > L. The only change is that ( [4.80 ) is now replaced by 



($ L , a (x)$ Lib (x / )) t = W L)a6 (|x-x'|;0 = 6 ab W L (\x-x'\;t), 



(5.118) 
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where Wl = J J IWI is now cut off at distance scaie L. Wl{0] t), its derivatives, and all 
relevant quantities constructed from Wl are ultraviolet finite. The distribution of zeroes, 
or line zeroes, of $l is given in terms of Wl as in the previous section. 

Choosing L = M~ l solves all our problems simultaneously. At wavelengths k~ l < L 
(i.e.k > M) the field fluctuations are oscillatory, with time scales 0{M~ 1 ). Only those 
long wavelength fluctuations with k^ 1 > L have the steady exponential growth that can 
lead to the field migrating on larger scales to its groundstates. Further, as the field 
settles to its groundstates, the typical vortex thickness is 0(M _1 ), and we only wish to 
attribute one zero to each vortex cross-section (or one zero to each monopole) . By taking 
L = M _1 we are choosing not to count zeroes within a string, apart from the central 
core. We shall see later that it is not crucial to take L = M~ l exactly, but sufficient to 
take L = C^M" 1 ). 

We are now in a position to evaluate Pt[3>], or rather Wl(t; t), which we now write as 
WM{f]t) for t > 0, where 



1 + A(k)(cosh(2tt(k)t) - 1) 



(5.119) 



1 + I,4(fc)(e n «* - e - n «*) 2 
2 



(5.120) 



W M (r;t) = / 4 u ke lk ~C(k) 

J\k\<M 

or,equivalently 

W M (r;t) = [ 4 D ke^C(k) 

J\k\<M 

and calculate the density of zeroes accordingly. Even though the approximation is only 
valid for small times, there is a regime Mt > 1, for small couplings, in which t is large 
enough for e Mt ^> 1 and yet Mt is still smaller than Mt sp = 0(ln(l/A)) when the 
fluctuations begin to reach the spinodal point on their way to the ground-state manifold. 
In this regime the exponential term in the integrand dominates and 

W M (r-t) ~ -/ d D kC(k)A(k)e ik -* e m{k)t (5.121) 

2 J\k\<M 

= TTFTTTn f dtt dk k ' 1 e 2n{k)t C(k)A(k)e^ (5.122) 

2\2lX) u J|k|<Af 

For the moment we assume that M and mo are comparable. In our context of initial 
thermal equilibrium at temperature To this corresponds to beginning the quench from well 
above the transition if, as earlier, we identify mo with the thermal mass at temperature 
T , most simply approximated to one loop by 

m 2 = -M 2 (l-^j (5.123) 



where T c = 0(M/\/\) in the same approximation. The definition ( |5.123| ) needs some 



qualifications, but for the moment we accept it as it stands . We shall consider the effect 
of modifying it later. For Tq > T c , and weak coupling C(k) of ( |5.116| ) is approximately 



C(k) = ^ (5.124) 
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The integral at time t is then dominated by the peak in k D 1 e 2n ^ t at k around ko, 
where 

D — 1 / / 1 



tk 2 = — 1 + 



2 -v--(m))- < 5 - 125 ' 

To check the consistency of the assumption that C(k) and A(k) are slowly varying in 
comparison to this peak we observe that, at the largest times of interest to us, t = 0(t sp ) 

k 2 = o(—)^M 2 (5.126) 

and thereby, for a quench from well above the transition, equally less than Mq. C(k) of 
( |5.124| ) can be approximated by 

C(k) = % (5.127) 

The effect of changing the initial thermal conditions is only visible in the 0(1/Mt) term. 
Since the overall normalisation of W(r; t) is irrelevant, both C(k) and A(k) m 1 can be 
factored out. For weak coupling we recover what would have been our first naive guess 
for the coarse-grained correlation function ((j)L,a( r ,t)4>L,b( r \t)) based on the growth of 
the unstable modes o (k, t) ~ e n ( fc )* alone, 

W M (r;t)<x[ d D k^e m{k) \ (5.128) 

J\k\<M 

provided we have not quenched from too close to the transition and t > M~ x (preferably 
t ^> M" 1 ). It is in this sense that our conclusions are insensitive to the initial conditions. 
However, we might hope that this insensitivity extends to non-thermal initial conditions 
as long as we are far from the transition. 

We note that if we had coarse-grained to longer wavelengths |k| < M < M then, once 
Mt > M 2 /Mq the dominant peak lies inside the integration region and the results are 
insensitive to the value of the cutoff. Provided M/Mq ~ 1 this will be the case because 
of ( [51261) . 

On approximating the peak around ko by a Gaussian (which assumes this insensi- 
tivity to coarse-graining on somewhat larger scales) we see that, for monopoles in two 
dimensions Wm{t", t) has the form of ( |4.89| ). 

W M (r;t) ex J dkJo(kr) exp{-^(k - k (t)) 2 /(Ak (t)) 2 } } (5.129) 
where ko(t) is given by ( |5.125| ), and 

Afco(t) 1 



k)(t) y/2' 



(5.130) 



Further, for the more interesting case of vortices in D = 3 dimensions, Wm(t", t) can 
be approximated by ([4.103 ), 

W M (r;t) oc J dksincikr) exp{-^(k - k {t)) 2 /{Ak {t)) 2 }, (5.131) 



25 



where 



is also large. We interpret (|5.130| ) and ( |5.132|) as implying large variance in the separation 



of the zeroes that are defining our defects. 

Although the forms (|5.129|) and (|5.131|) represent the spread in wavenumber well they 



do so at the expense of a correct description of the long wavelength modes. WM(r; t) of 
( |5.128| ) can be compared directly to the empirical form of (|4.104|) by writing it (up to 



the usual arbitrary normalisation) as 

W M (r;t)~ I dksmc{kr)(-^—] 2 e 2n ^ 1 . (5.133) 

J\k\<M \k (t)J 



In the terminology of ( |4.104j ) this shows that, for long wavelengths, the power is n = 0, 
determined entirely by the radial k ' 1 behaviour in D = 3 dimensions. Insofar as 
distributions of vortices are determined by the power in long wavelengths the peak at 
k = ko(t) may, to some extent, be approximated by the cutoff at ko(t) of (|4.104|) . We 
shall return to this later. 



Rather than evaluate correlation functions from (|5.128|) directly it is convenient to 



approximate it differently. It is not difficult to see that, for the values of Ak (t)/ko(t) 
given above, WM{r;t) approximately falls monotonically in r to zero from above (any 
oscillatory behaviour is of very small amplitude). It is a good approximation if, in ( |5.128|) 
we expand fl(k)t as 

Q(k)t^Mt-^ (5.134) 

from which 

W M (r; t) ~ e 2Mt [ d D k e ik x e- kH/M . (5.135) 

J\k\<M 

This approximation maintains the peak in the integrand at ko(t) of ( |5.125 ). Taking it as 
it stands gives 

f(r; t) = exp{-r 2 M/4t}, (5.136) 

independent of dimension, on dropping the upper integration boundQ. As such it cor- 
rectly reproduces the small-r behaviour of f(r;t) (and f°(r;t) of (|4.106|) ) and, although 
the simple Gaussian behaviour breaks down for r 2 M/4t ^> 1, the rapid large-r falloff is 
qualitatively correct. 

We conclude with one more observation. Even in the disordered state for t < 
there were zeroes (or lines of zeroes) induced by thermal fluctuations. Effectively, the 
equal-time two-point correlation function for t < is 

W M (r;t) ~ T / jPk f "* 2 = / d D k e ik ^C(k). (5.137) 
j|k|<M k z + rriQ J\k\<M 



3 This is not a question of dropping coarse-graining but merely approximating the integral. 
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Because of equilibrium there is no time- dependence in WAf(r;t) of ( |5.137| ). There is an 
implicit cutoff in ( |5.137|) at the thermal wavelength (3q = Tq , but we are not interested 
in zeroes on a smaller scale than the same field correlation length M -1 that determines 
defect size. In fact, assuming itlq m M, it makes little difference whether we coarsegrain 
to |k| < M or to |k| < mo- 

Whatever, the surface density of zeroes is 

n{t) = 0{ml) = 0{M 2 ). (5.138) 

However, we should not think of these zeroes as the precursors of the defects that appear 
after the transition. They are totally transient. We see this by the presence of oscillatory 
factors e tul ^ At in the two-point correlation function at unequal times t,t + At. In fact, 
the term ( 5.137|) is the first term in (|5.119 ) and (15.120 ), in fact the only term in either 
of these expressions when t = 0. The reason for rewriting (|5.119| ) as ( |5.120| ) was to 
discriminate between this term (the 1 in the square bracket of ( |5.119|) ) and the term 
—A(k) ~ — 1, in the same bracket, whose origin is the interference between exponentially 
increasing and decreasing terms in the mode evlution and has nothing to do with thermal 
fluctuations directly. Since the normalisation of WM{f)t) has no effect on the density 
and distributions of zeroes and line zeroes these thermal fluctuations are suppressed, 
relative to the long wavelength peak that shows the growth of domains, by a factor 
0(e~ 2Mt ). Thus, although the thermal fluctuations remain as strong in absolute terms, 
their contributions to the counting of zeroes, and hence the creation of defects, vanishes 
rapidly, which is why we did not include them earlier. 

Similarly, had we chosen to coarse-grain on a somewhat shorter scale |k| < Mo, where 
Mo > M, then W(r; t) would have acquired some (finite) oscillatory terms from ( |5.114 ) 
that in turn would have been suppressed, relative to the long wavelength peak, by a 
factor 0(e~ 2Mt ). For large enough times t sp they would play no contribution, provided 
M/Mq ~ 1. Taken with our earlier observations on coarse-graining at larger scales than 
M -1 we see that coarse-graining on a scale comparable to the thickness of cold defects 
is all that is required, with fine-tuning being unnecessary. 



6 Slowing the Quench 

Our adoption of an instantaneous quench, in which m 2 (t) changes as in (|5.112p is obvi- 
ously unrealistic. Any change in the environment requires some time r to implement. In 
the context of the flat spacetime free-field approximation that we have adopted so far 
we should replace m 2 {t) of fl5.112j ) by an effective {mass) 2 that interpolates between m\ 
and its cold value — M 2 over r. If r were so small that Mr < 1 we expect no significant 
change, since the field would not have responded in the time available anyway. However, 
once 

1 < Mr < Mt sp (6.139) 

there will be an effect. As the wavenumber k of the field modes increases towards M the 
time available for their growth is progressively reduced. However, the long wavelength 



27 



modes have the same time to grow as before. As a result, the peaking of the integral of 
W(r; t) will occur at a smaller value of k for the same time t after the quench is begun, 
leading to a lower density of defects than would have occurred otherwise. Alternatively, 
the field is correlated over larger distances than would have happened otherwise. 

To see how this occurs quantitatively it is not necessary to go beyond simple approx- 
imations, in the absence of any compelling reason to make a specific choice for m 2 (t). 
[Such a reason occurs if we have a definite spacetime metric driving the transition, as 



happens in inflationary models [13 
Consider the behaviour 



mr(t) = 777, , t < 



= — M — , < t < t 

T 

= -M 2 , r<t (6.140) 

in which the quench is begun at time t = and completed at time t = r. The rate of 
the quench, r _1 is assumed large, r < t sp . That is, 

1 < Mr < ln(l/A), (6.141) 

so that the quench is complete before the field has experienced the turnup of the potential 
towards its minima. As r — > we recover the instantaneous behaviour examined earlier. 
If the behaviour presented in (|6.140| ) seems a little artificial in how the behaviour 



before t = and after t = are connected we note that, from our previous discussion, 
the behaviour of m 2 (t) for t < is largely irrelevant, as long as m 2 is comparable to M 2 . 
Only the behaviour for t > will be important. 

The calculation of f(r; t) now requires a proper solution of the mode equation ( |5.111| ) 



for m 2 (t) of (16.140]) ) in terms of which W(r;t) is again given by (|5.113|) . This is not 



possible to derive analytically, but it is not difficult to make analytic approximations. 
As before we coarse-grain W(r;t) to eliminate oscillatory modes, bounding |k| by M. 
Whereas modes of all wavelengths |k| < M begin to grow exponentially at the same time 
t = for an instantaneous quench (r = 0), for r ^ modes of wavenumber k do not 
begin to grow until time t T (k), at which 

k 2 = -m 2 (t T (k)) (6.142) 

For the choice of m 2 above in (|6.140| ) 



Tk 2 

tr{k) = W2 (6.143) 

As we noted earlier, although the very long wavelength growth is activated at t ~ 0, the 
shorter the wavelength (but still |k| < M) the shorter the time the modes have available 
before their switchoff at t = t sp = O^M' 1 ln(l/A)). 
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A crude, but helpful approximation, for times Mr < Mt < ln(l/X) is to mimic this 
effect by modifying ( |5.128| ) as 

W Mr (r;t)~ [ 4D k ^ e 2n(k)(t-t T (k)) : ( 6 .i44) 

J\k\<M 

In the further approximation of (|5.134 ), in which we expand in powers of k, 



W M Jr;t) ~ / 4D ke i^ e (t-t T (k))(2M-k*/M) 

J\k\<M 

[ d D k^e- k ^ t+2T ^ ML ' 

J\k\<M 



J2Mt 



l + 0( 



M 3 ' 



(6.145) 



On comparing WM, T (r;t) of ( |6.145| ) to that of ( f5.135| ) it follows that, at the level 
of the Gaussian approximation, the effect of slowing the quench is to replace f(r;t) of 
(Vm) by 

f T (r; t) « exp{-r 2 M/4t r } (6.146) 

where 

*^ t + 2T+0 {w^j} (fU47) 

That is, in the first instance the effect of slowing the quench is no more than to reproduce 
the behaviour of an instantaneous quench, displaced in time by an interval 2r. A more 
careful calculation would give slightly different answers, but the qualitative conclusion 
that slowing the quench reduces the defect density and hence makes the system look as if 
it had begun earler, is correct. In the absence of any reason to make a particular choice 
of quench, the result ( |6.147| ) is adequate for our purposes. 

7 Freezing in the Vortices 

We understand the dominance of wavevectors about k\ = M/t T in the integrand of 
W(r;t) as defining a length scale 



Ut)=0(^U/M), (7.148) 
once Mt > 1, over which the independently varying fields <ft a are correlated in magnitude. 



To be specific, let us take £ T (t) = 2yt T /M, whence 

f T (r;t) = exp{-r 2 M/4t T } 

= exp{-r 2 /e T (t)}. (7.149) 

With this definition, the number density of line zeroes at early times is calculable from 
( g]flj ) as □ 



11 



The density of monopoles in D=2 dimensions is identical. 
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both for monopoles and vortices, permitting us to interpret £ T {t) as a correlation length 
or, equivalently, a domain size for phases. Although the zeroes and line zeroes have yet 
to freeze in as defects, this density of one potential defect per few correlation areas is 
commensurate with the Kibble mechanism cited earlier. 

The zeroes that we have been tracking so far cannot yet be identified with the vortices 
(and monopoles) which provide the semiclassical network on the completion of the tran- 
sition because the field has not achieved its groundstate values. In fact, and of greater 
importance at this stage, it is not even uniform in magnitude. Insofar as a classical 
picture is valid, the thermal fluctuations have kicked the field off the top of the upturned 
potential hill in different directions at slightly different times. The work of Guth and 
Pi||13|| shows that the variance in this time is At = 0(M~ 1 ). Thus, even if the poten- 
tial were a pure upturned parabola and the fields were to stop growing instantaneously 
at the moment that they reached the spinodal value V"(<ft sp ) = (defined in terms of 
the non-Gaussian physical potential) it would take a further time O(At) before the field 
caught up in all places. 

An instantaneous halt to an increasing field growth is a huge oversimplification. More 
realistically, as the fields approach their spinodal values the effective (mass) 2 driving the 
expansion of the unstable modes decreases and the expansion of the domains slows. The 
energy of the fields in the long wavelength modes will cause an overshoot towards the 
potential bottom that, in the absence of dissipation, will be followed by a rebound, and 
perhaps a repeat (or many repeats) of the whole cycle[|16|. The details will depend on 
the dissipation that the long wavelength modes endure. This dissipation is necessary to 
enable the fields to relax to their ground-state values, and for realistic fields has several 
sources. In particular, the act of coarse-graining induces dissipation HT7|, as would the 
presence of other fields [H| Q. 

In this excursionary paper we shall just consider the rudiments of the effect of the 
field self-interaction on slowing down and stopping domain growth in an approximation in 
which dissipation is assumed to set in rapidly at the spinodal field values. At the simplest 
level this is, in some respects, like that of slowing the quench that we discussed previously, 
except that it occurs at the end of the period of interest, rather than the beginning. As 
there, the long wavelength modes are the least affected, the shorter wavelength modes 
(but still with k < M) having less time to grow. For this reason the density of zeroes, 
now identifiable as defects, is reduced. 

As we have said, a full analysis is beyond the scope of this paper but, as with the 
slow quench, it is possible to provide a primitive approximation in which we can see 
this explicitly. This is enough for our present purposes, for which the power in the long 
wavelength modes will turn out to be the most relevant property that we need, and we 
hope that the crudity of our approximations leaves this untouched. 

If we wish to retain the Gaussian approximation for the field correlation functions the 

12 A further source of dissipation is an expanding metric, inevitable in cosmological models. This is 
beyond us here. 
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best we can do is a mean-field approximation, or something similar^. In this approxi- 
mation equations of motion are linearised by the substitution 



(7.151) 



and similarly for <j) 2 . Because of the diagonal nature of W a & this means that 



% + ft) 



(7.152) 



W(r; t) still has the mode decomposition of ( |5.113| ), but the modes U^ k now satisfy the 
equation (ignoring subtractions) 



^ + k 2 + m 2 (t) + 4A J 4 D pC(p)U+ p (t)U- p (t) 



(7.153) 



A detailed discussion of (|7.153|) will be given elsewhere, but a rough estimate of the 
effects of the interactions can be obtained by ignoring self-consistency and retaining only 
the unstable modes in the integral. For simplicity we revert to an instantaneous quench 
m 2 (t) = —M 2 (t) for t > 0. On using the definition for A in Q5.109 ), we might approximate 
( [7.1531 ) in turn by 

il + k 2_ M 2 + M 2 e 2M(t- W ) 

at 2 



0. 



(7.154) 



where, further, we have ignored all but the dominant exponential behaviour due to a 
free-field roll. Note that the A-coupling strength has been absorbed in the definition of 
t sp , most justifiable for small A. Equation ( [7.1541 ) has the correct qualitative behaviour 
in that there are growing modes only for t < t sp , after which all modes are oscillatory. 
In this approximation a mode U^k can only grow until time tf(k), defined by 



£2 = M 2 (l- e 2M{tf{k) ~ tap) 



(7.155) 



As anticipated, shorter wavelengths have less time to grow. Linearising ( 7.155| ) in the 
vicinity of t sp as 



k = 2M-\t f (k)-t 



sp I 



(7.156) 



shows that, beginning from an instantaneous quench, the mode with wavenumber k stops 
growing at time 

before the modes of longest wavelength stop. 

Instead of the instantaneous quench, now let us reintroduce a slower quench that takes 
time r to implement, along the lines of ( p. 140 ) and the approximations that followed from 



13 But not a large-iV 0(N) limit, since we only have simple defects in D — N and D — N 
dimensions. 
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it. At the simplest level, the effect of the back- reaction on the field correlation function 
W M , T (r; t) of (|6.144|) is to replace it by 



where 



J\k\<M 



t*{k) « t T {k) + At{k) 
rk 2 k 2 



+ 7777?- (7-159) 



M 2 2M 3 ' 

In the further approximation of (|5.134 ), in which we expand in powers of k, 



W M Jr,t)~ / ^ jfce *.x e (t-r(fc))(2M-fcVM) i (7160) 

J\k\<M 

From this we derive an f(r; t) of the form 

/(r; *) « exp{-r 2 M/4t T } (7.161) 

where t T is now given by 

t T ^t + 2r + ^. (7.162) 

rather than by 147| ) . We stress that, even at this level of approximation (in which any 
term can be multiplied by a coefficient 0(1)), we have assumed sufficiently weak coupling 
and a rapid enough quench that 

r + — <t SJ> (7.163) 
2M p K ' 

so that even the shortest wavelength considered has some time to grow. The details of 

( [7.1611 ) and ( [7.162] ) should not be taken too seriously, but the qualitative result that the 

effect of slowing the quench is to give a time-delay of 0(t), and the effect of back- reaction 

is to give a time-delay of 0(M _1 ) is as we would expect. A more realistic calculation |i~6f 

based on (|7.153|) would show that the field overshoots its spinodal values, retreats and 



overshoots again in damped oscillations. However, the small-r behaviour of ( |7.161| ) is 
good enough for our immediate purposes. We shall return to the full equation ( |7.153|) 
later. 



8 Densities and Correlations from a Quench 

In the approximation given above the effect of slowing the quench (and taking the back- 
reaction into account) is just to increase the correlation length and thereby reduce the 
density of zeroes to one appropriate to an instantaneous transition that happened earlier 
and in which the field spontaneously stopped growing. The only way the zero-density can 
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decrease without the background space-time expanding is by zero-antizero annihilation. 
For monopoles this is simple removal of zeroes by superposition. For vortices it will 
include the collapse of small loops of zeroes. The end result is to preserve roughly one 
string zero per coherence long held belief for whatever mechanism. 

Let us perform all our calculations at t = t sp , for which £ T (t sp ) = £ sp . For small 
enough A, the density at t — t sp is, with our previous caveats about coefficients, 

1 1 1 M 2 

n(t sv ) = -4- ~ —7 7 < m2 , ( 8 - 164 ) 

1 sp) 7T e sp 4tt (Mt sp + 2Mt + 1) y J 

We have retained the last term in the denominator as a reminder that the approximation 
may work even when Mt sp > 1, but not too large. This shows that the line zeroes (or 
monopole zeroes) only create a small fraction of space as false vacuum. Equivalently, 
the typical separation of line zeroes (or monopole zeroes) is significantly larger than the 
thickness of a cold vortex (or monopole) once they have frozen in. This is one reason 
why we can begin to consider these line zeroes (and point zeroes) as serious candidates 
for vortices (and monopoles) even though the field has not fully relaxed to its ground 
states. 

A further reason is that, even though the spinodal values of field are at a fraction of 
1 / V3 of the distance to the bottom of the potential, this distance is very much larger 
than the initial field fluctuations. In consequence, the fluctuations in the field are now 
too small to create new lines and points of false vacuum of any substance and we can 
safely say that we have defects. To see this, we observe that, initially, the probability 
Pt[<E>], now independent of t, is expressable as 

p [$\=Ne- fhHm (8.165) 

N is a normalisation factor and H permits the expansion in f3o m o [0 m terms of Sq of 

mm) 

1 r ft) S 1 \ 2 
H[$} = S m - Y / J d D x + O(0X). (8.166) 

With calculational simplicity in mind we restrict ourselves to high initial temperatures 
(/?o^o ^ 1) f° r which it is sufficient to retain only the first term, and take Ao = 0. 

As with zeroes, fluctuations are defined with respect to some length scale L, say. The 
fields $>L,a coarse-grained to this length are defined in ( [5.117 ). The probability Pl{&) 
that $l = $ can be written as 

p L ($) = N exp(-$ 2 /2(A L $) 2 ) (8.167) 

where (A^$) 2 ) is the coarse-grained two-point function 



(A L $) 2 = W L (r;0)~T [ fh 



!<£,-! k 2 + m 



d D kC(k) (8.168) 

Ik^L- 1 
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for C(k) of ( |5.116| ), compatible with ( |5.119| ). For mo ~ M it follows |I5] that, for L = 
0(M~ l ) in D = 3 dimensions 

(A L $) 2 ~ AMT (8.169) 

where A ~ 10 _1 . 

The condition that there be no overhang is thus 



(A L $) 2 « <P 2 



(8.170) 



There is no difficulty in satisfying flS.170| ) for small coupling, for which it just becomes 
A< 1. 

Although the strings of zeroes (or monopoles) are far apart they still do not yet quite 
provide the semiclassical network of defects that can be used as an input for numerical 
simulations since the fields have to relax dissipatively from their spinodal values to the 
true minima of the potential along the lines suggested earlier. Because they are too 
costly in energy to be produced by fluctuations, there will be some defect annihilation, 
but no creation, in this final phase of freezing in the defects and the n(t sp ) calculated 
previously is an overestimate of the string density at the end of the transition. However, 
if it were the case for vortices that all the string was in loops it is difficult to see how 
infinite string could be created if all that happens is that string is removed from the 
system, although it cannot be precluded. Since that is the main question that we shall 
be addressing here we adopt the simplest assumption that, even if the density n(t sp ) is 
an overestimate the distribution of strings (i.e. the fraction of strings in loops, the index 
for length distributions) is approximately unchanged by the final freezing in. That is, 
the one-scale scaling regime remains valid till freeze-in. 

In that case the distribution of strings will then be as above for £ = £(t sp ), while 



A£/£ of ( |8.183 ) remains (albeit approximately) the variation in domain size over which 
field phase is correlated. Direct attempts to determine distributions from density corre- 
lation functions are difficult. For example, the simple analytic form of (|7. 149|) enables 
us to calculate the density correlation function C(r; t sp ) for monopoles in the plane, and 
A(r;t sp ) and B(r;t sp ) for vortices in three dimensions, up to exponentially small terms 
in Mt T (t sp ), as 



C(r) = A{r) 



g 2r / £ sp 



and 



7T 2 e p (l_ e -2rVC? p )2 

2 e 



><sr> 



<0, 



Br) 



n 2g p (I - e -2rV& p . 



>0, 



(8.171) 



(8.172) 



and we have dropped unneccessary time labels. 

In units of n 2 (t sp ), the anticorrelation is large. In particular, a calculation of C(r) or 
A(r) for defect separation r gives 

,2 



A{r) = n\t sp ) 



1 + 

3 SP 



(8.173) 
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when r < £ sp , as happens for the single mode case. Again, as in the single-mode case, 



B(r) = n 2 (t sp ) 



|f 1+0(1) 



> 0. (8.174) 



Although ( |8.172[) is not wholly reliable for r £, the suggestion that C or A and B 
fall off very fast is qualitatively correct, showing that £ sp indeed sets the scale at which 
strings see one another. The best we can say is that, taken together these suggest a 
significant amount of string in small loops. 

We have more success in using the correlation functions to determine the variance 
in vortex winding number N$ through an open surface S in the 1-2 plane, which we 
take to be a disc of radius R. There is a complication in that, as can be seen from 



( p.32j ), P3 counts zeroes weighted by the cosine of the angle with which they pierce S, 
whereas winding number counts zeroes in S without weighting. The problem is therefore 
essentially a two-dimensional problem. Using the results of ( 2.28|) onwards, it is not 



difficult to see that, given the short range of the correlation functions, for the density of 
zeroes 

(AiV)| = 0(n 2 e sp R) = 0{R/i sp ). (8.175) 

Since (AN)s is (Aa)s/27r, where (Aa)s is the variance in the field phase around the 
perimeter dS of S, this means in turn that 

(Aa) 2 s = 0(R/£(t sp )). (8.176) 

If this were equally true for nonrelativistic vortex systems, then (Aa)| measures the 
variance in the supercurrent produced by the quench [0] and such behaviour is measurable, 
in principle. In practice we do not yet know how to construct the coefficient of R for 
non-relativistic fields, but see Refs.|2(| |21|| . 

Returning to the main problem of the length distributions of vortices in three di- 
mensions, it has not yet proved possible to turn expressions for the Cy (i.e. A and B) 
directly into statements about self- avoidance, fractal dimension, or whatever is required 
to understand the resulting string network. Fortunately, as a temporary expedient, we 
can work indirectly by adapting the numerical results of p|, based on Gaussian fields 
with the two-point function W n (r;t) of ( f4.104| ) that we considered earlier. For n = 



this has the same long- wavelength behaviour as our dynamical W(r; t) of Q5.133 ) and 



( |5.135| ). The results of || for n = reproduce those of || on a cubic lattice, with its 
preponderence of open string. This leads us to believe that the simple dynamical model 
that we have proposed above (given its assumptions for the freezing in of domains) would 
also lead to a large fraction of open string. 

Extending the simulation to general n for W n (r;t) of (|4.104p showsH that, in order 



to decrease the amount of string in open string, it is necessary to decrease n. However, 
although the fraction of infinite string does diminish as n becomes negative, infinite 
string only seems to vanish in the pathological limit n — > —3, at which there is an 
infrared divergence. We see from (|4.105|) that the density of string vanishes in the same 
limit. As long as there is a nonzero density of string, then some of it is infinite. 
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The only way to increase the power in long wavelengths in our dynamical model is to 
change C(k) by imposing different initial conditions. Altering C(k) from C(k) = 0(k°) to 
C(k) = 0(k n ) leads to a change in the power of the long wavelength modes from n = 
to n, corresponding to a modification of WM^r; t) of ( |5.133|) to 



W M {r;t)(x dk sinc(kr) k 2+n e 2Q(k)t . (8.177) 

J|k|<M 

The effect of changing n can be seen qualititively by expanding exp{Q(k)t} as in 
( |5.134| ). WM{r;t) of (18.1771 ) can then be evaluated explicitly, giving f {n) {r-t) as the 



confluent hypergeometric function 

/W(r;t) w expi-r'M/At}^—;^; — 
/n + 3 3 r 2 M\ 

on dropping the upper bound in the integral. 

The densities obtained from the more complicated (|8.178|) are exactly those of the 
more simple W n (r; t) given in ( |4.105| ), if we identify k 2 as M/2t. Further, for n > —2 the 
integrand of W(r; t) in ( p.!77| ) shows a peak whose width increases as 



A&oc^==, (8.179) 

implying a variation in domain size that also increases asn-* —2, just as for the simpler 
case. 

Since i-Fi(0; 3/2; r 2 M/4t) = 1 identically, n = is a special case in which we reproduce 
the Gaussian of ( |5.136| ). Although this seems at variance with the behaviour given in 
( |4.108j ) for the simply cut-off power distribution, the exponential behaviour was never 
supposed to be valid for very large r. Despite that, for n < the similarity is good, 
with the same power-law falloff for n — — 1 and n = —2, leading to anticorrelations over 
increasingly longer ranges. 

Thus, if we can force C{k) to have negative n, the numerical simulations of suggest 
that there will be more string in loops, and less in infinite string. The reason why we 
have n = white noise in ( |5.133| ) is that, when quenching from well above the transition, 
C(k) w Tq/ttLq is constant for the dominant integration region. In attempting to alter C(k) 
we expect that, if we were to begin closer to the transition where the initial fluctuations 
are larger, we can make n negative. On beginning a quench from thermal equilibrium, the 
most extreme case is one in which we start from a temperature so close to the transition 
that the effective mass ttlq of (|5.112 ) for t < is approximately zero. Let us first suppose 



that we could set it identically zero. Then 



C(k) = f 2 . (8.180) 
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rather than ( 5.127|) . The k~ 2 behaviour cancels the k 2 behaviour coming from the radial 
momentum integration and the power in long wavelengths is increased from n = to 
n = —2 in ( |S.177| ). We believe that this is the most negative that n can become. Even if 
we had performed a more honest calculation of (|7.153|) there is no way for the interactions 
to introduce the singular infrared behaviour of n = —3 that is necessary to prohibit the 
production of strings. The simulations of || suggest an approximate halving of the 
fraction of open string for n = —2 from its original 80% (on a cubic lattice). These 
numbers are not to be taken too seriously since, as we commented earlier, they depend 
on the type of lattice^ (although the presence of open string does not). 

As we had noted above, for such a C(k) there is no peaking in the integrand around 
any wavenumber, and so no length that characterises a domain size, although the string 
density is nonzero. From ( |8.168| ) we see that the same C(k) of (|8.180 ) also determines the 
initial field fluctuations. However, the effect there is small, contrary to naive expectation, 
provided we coarse-grain to the scale appropriate to defects, k < M. This is the relevant 
scale, rather than the thermal wavelength (3 Q . On coarse-graining to k < M the field 
fluctuations, given by ( |8.168| ), still satisfy (A L $) 2 ~ AMT , and the initial fluctuations 
will not populate the ground states for weak coupling. 

Our ability to make n negative survives more realistic initial states. The situation is 



more complicated in that the Gaussian approximation based on ( 5.110 ) breaks down|15| 
at temperatures T closer to T c than 



rp2 

~ J>2 



0(A). 



(8.181) 



Thus the minimum value of at which our approximations have any hope is 



m 



0(XM ). However, this is small enough for flB.18U| ) to be a good approximation, provided 



A ln(l/A) < 1, 



(8.182) 



as follows on implementing the approximation (|5. 134 ) . In practice quenching from so 
close is artificial, and we should consider n = —2 as an unattainable lower bound. 

However, this is not the end of the story. The regular lattices of and are 
artefacts of the calculational scheme and may, of themselves, predict more infinite string 
than is present if they try to emulate a continuous transition. As we saw earlier, we do 
not have a regular domain structure in our model but have domains with a large variance 



Ak 



0(1) 



(8.183) 



independent of time. Whatever the details, the dispersion in £ is large. 

Although a regular lattice does not imply exactly regular domains (which we under- 
stand as characterising the spaces between the defects) for high densities of the kind we 
have here they do imply a greater regularity than we have. In numerical simulations of 
string networks the inclusion of variance in the 'lattice' cell size showsH that, the greater 



37 



the variance, the more string is in small loops. This can be understood in the following 
way. The strings generated by phase separation on regular lattices are known to behave 



like random walks in D = 3 dimensions to a very good approximation [22]. On a rectan- 
gular lattice the fraction of string in loops is only about 20% ||. Increasing the variance 
of the domain size increases the 'target area' that a string must hit for a section of string 
to be deemed closed. The probability of finding loops therefore increases. The variance 
of (|8.183|) , if it could be carried over to || as it stands, suggests much more string (e.g. 



twice or more) in small loops. However, it is not easy to marry the somewhat different 
distributions of this simulation to ours since ours is one of overlapping domains (because 
of the continuous transition) rather than a domain 'bubble' picture more appropriate to 
first-order transitions that is more easily accomodated by ||. 

The end result is that there always seems to be infinite string from a continuous 
transition, even if less than we thought, although its confirmation will require numerical 
analysis along the lines of M. This will be performed elsewhere when we have a better 
understanding of how back-reaction stops domain growth [p^] than the simple approxi- 
mation presented here. 



9 Conclusions 

In this paper we have shown how global 0(2) vortices (and monopoles) appear in a 
relativistic theory, at a quench from the ordered to disordered state, as a consequence of 
the growth of unstable Gaussian long wavelength fluctuations. 

Assuming a simple freezing of defects the resulting string (and monopole) configura- 
tions scale as a function of the correlation length £ = 0(\/tP), for some effective time t T , 
at a fraction of a defect/correlation area. The label r characterises the time it takes the 
defects to freeze in. The slower the quench (for weak coupling), the larger t T , and hence 
the lower the defect density. For quenches from well above the transition the results are 
insensitive to initial conditions, which have been assumed here to be thermal equilibrium, 
but possibly the result is more general. This is compatible with the Kibble mechanism 
for vortex production due to domain formation upon phase separation by white noise 
fluctuations (and similarly for monopole production). For a weak coupling theory the 
domain cross-sections are significantly larger than a vortex (or monopole) cross-section 
at the largest times for which the approximations are valid. Moreover, there is a large 
variance in their size £, with A£/£ = 0(1). Numerical simulations that relate loop dis- 
tribution to the variance in domain size of the long wavelength modes suggest that there 
is likely to be more string in small loops than we might have anticipated on the basis of 
simulations on regular lattices, although the match with the simulations is not exact. 

As we quench from closer to the transition and effectively reduce n we expect that 
the fraction of string in loops increases (as the density of string decreases), but even then 
it is never possible to enhance long wavelength fluctuations to an extent that there is no 
infinite string after a continuous transition. 

We stress the crudity of some of the approximations that we have made. There 
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is little difficulty, in principle and in practice, in doing somewhat better. However, 
without any specific choice of initial conditions and mass evolution m(t) to guide us, and 
without better-matching numerical simulations, generic results at this level are sufficient. 
Improvements, motivated by particular cosmological models, are being considered. 
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